Radion Stabilization by Stringy Effects in General Relativity 



Subodh P. Patil 1 ) and Robert Brandenberger 2 ' 1 
1) Department of Physics, Brown University, 
Providence, RI 02912, USA 
2) Physics Department, 
McGill University, 3600 rue Universite, 
Montreal, QC H3A 2T8, CANADA 
E-mail: patil, rhb@het.brown.edu 

We consider the effects of a gas of closed strings (treated quantum mechanically) on a background 
where one dimension is compactified on a circle. After we address the effects of a time dependent 
background on aspects of the string spectrum that concern us, we derive the energy-momentum 
tensor for a string gas and investigate the resulting space-time dynamics. We show that a variety 
of trajectories are possible for the radius of the compactified dimension, depending on the nature of 
the string gas, including a demonstration within the context of General Relativity (i.e. without a 
dilaton) of a solution where the radius of the extra dimension oscillates about the self-dual radius, 
without invoking matter that violates the various energy conditions. In particular, we find that in the 
case where the string gas is in thermal equilibrium, the radius of the compactified dimension dynam- 
ically stabilizes at the self-dual radius, after which a period of usual Friedmann-Robertson- Walker 
cosmology of the three uncompactified dimensions can set in. We show that our radion stabilization 
mechanism requires a stringy realization of inflation as scalar field driven inflation invalidates our 
mechanism. We also show that our stabilization mechanism is consistent with observational bounds. 



I. INTRODUCTION 

In the early days of string cosmology, it was realized 
that superstrings had an effect on space-time dynamics 
that was qualitatively quite different from that of parti- 
cles or fields. In particular, it was realized that string 
winding modes could provide a confining mechanism for 
certain compact directions in such a way as to allow only 
three spatial dimensions to grow large Q. Key to this 
realization are the T-duality of the spectrum of string 
states, and the fact that the background is described by 
Dilaton Gravity, and not by General Relativity with a 
fixed dilaton (this is crucial in order that the background 
equations obey the T-duality symmetry). The arguments 
of were put on a firmer basis by the analysis of 0] (see 
also 0). 

Starting point of the considerations of |lj is the as- 
sumption that all spatial dimensions begin at close to 
the self-dual radius (the string scale), and that matter 
consists of a hot gas of string states. The considerations 
of p] were more recently applied to "brane gas cosmol- 
ogy" HI j a scenario in which the initial string gas is 
generalized to be a gas of all brane modes. It was shown 
that given the hot dense initial conditions assumed in , 
the string winding modes are the last modes to fall out 
of equilibrium and thus dominate the late time dynam- 
ics. Hence 0, the inclusion of brane degrees of freedom 
does not change the prediction that only three dimen- 
sions to grow large. The dynamical equations describing 
the growth of the three dimensions which can become 
large were solved in @ (see also @). In 0, it was shown 
that isotropy in these large dimensions is a consequence 
of the dynamics. In |8J, it was found that if both the mo- 
mentum and winding modes of the strings are included in 
the dynamical equations, the radius of the compactified 



dimensions is stabilized at the self-dual radius. More pre- 
cisely, the expansion of the three large dimensions leads 
to damped oscillations in the "radion" about the self-dual 
value. Thus, in the context of a background described by 
Dilaton Gravity, radion stabilization is a natural conse- 
quence of brane gas cosmology 28] . 

At the present time, however, the dilaton is most likely 
fixed (see, however, [llj for an alternate scenario). Thus, 
it is of interest to explore how the inclusion of string 
(and brane) winding and momentum modes influences 
the dynamical evolution of the radion in a background 
space-time described by General Relativity (GR). There 
is another motivation for studying this issue. Another 
corner of the M-theory moduli space is 11-d supergravity. 
In |12| it was found that a brane gas in this background 
also admits a region in the phase space of initial condi- 
tions in which only three spatial dimensions can become 
large, although this corner may not be consistent with 
holographic entropy bounds 01 ( see a ls° Q where the 
considerations in this corner of M-theory moduli space 
was extended to spaces with more general topologies). 
Motivated by these considerations, we in this paper study 
a simplified problem, namely the questions of how a gas 
of winding and momentum modes of strings winding one 
compactified spatial dimension (taken to be a circle) ef- 
fects the evolution of the radius (the radion). We start 
with initial conditions in which the three spatially non- 
compact dimensions are expanding. We find that the gas 
of string winding and momentum modes gives a natu- 
ral radion stabilization mechanism. Our approach is to 
consider the effect of strings on 5D space-time dynamics 
(with the extra spatial dimension compactified to a cir- 
cle) by adding the appropriate matter term to the stan- 
dard Einstein-Hilbert action. We will derive this term 
shortly (see also for a similar derivation). The result- 
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ing energy-momentum tensor leads to a novel behavior 
when inserted into the Einstein equations. We will find 
that we can generate a non-singular bouncing solution for 
the radius of the compactificd dimension in the context 
of GR (without a dilaton) while respecting the Domi- 
nant Energy Condition for the matter content. Specifi- 
cally, the radion performs damped oscillations about the 
self-dual radius. Initially, we study a pure state of mat- 
ter with specific quantum numbers obeying the T-duality 
symmetry. However, we will find that we can rather nat- 
urally extend the analysis to a gas of these strings in ther- 
mal equilibrium (with a bath of gravitons and photons), 
with the result that the radius of the compact dimension 
is dynamically stabilized at the self dual radius R = Va*, 
where 27ra' is the inverse of the string tension (see also 
|l6j | for a study of string gases in thermal equilibrium) . 

In addition, we find that our model evolves according 
to standard Friedmann-Robertson- Walker (FRW) cos- 
mology after the compact dimension has been stabilized, 
and that the resultant stabilization is incompatible with 
any subsequent inflationary epoch driven by a bulk scalar 
field (for string-specific ideas on how to generate inflation 
in brane gas cosmology see \Vi\). However this conclu- 
sion can be avoided if some form of stringy inflation is 
realised where strings are produced in re-heating. 

Before we can turn to any of this, we will have to ad- 
dress a question of principle concerning the string spec- 
trum in a cosmological context (this issue is also being 
studied in 18]). The question of formulating String The- 
ory in a time-dependent background is a current and 
active area of research. However, we are primarily in- 
terested in the behavior of strings in a background that 
evolves on a cosmological time scale. As can be seen from 
the FRW equations, the cosmological time scale H^ 1 
(where H is the Hubble expansion rate) is larger than 
the characteristic microscopic time <r -1 (where u 4 is the 
matter energy density) by a factor of m p i/a, where m p i 
is the Planck mass. Thus, away from singular epochs in 
the history of the Universe, the cosmological time scale is 
going to be many, many orders of magnitude longer than 
the characteristic time scale of the string dynamics, and 
hence we should be able to inherit many of the features 
of the string spectrum in a static space-time (with some 
obvious modifications). We justify this intuition more 
rigorously in the Appendix, but we feel that it might suf- 
fice at this point to remind the reader of the approximate 
irrelevance of a time dependent background for a much 
more familiar theory: Quantum Field Theory (QFT). 
Although quantum fields in curved spaces exhibit sev- 
eral qualitatively different features from quantum fields 
in flat spaces[23, we still manage to do a lot of sensible 
(and spectacularly successful) flat space-time QFT calcu- 
lations despite the persistent Hubble expansion of space- 
time. The reason for this is easy to see: the contributions 
to masses, to scattering amplitudes, to the structure of 
the Hilbert space of our theory, etc., that come from 
terms that depend on derivatives of the metric are in the 
present epoch highly suppressed and irrelevant. This is 



partly captured by the Adiabatic Theorem, which is the 
statement that given two systems with Hamiltonians that 
can be continuously interpolated, then in a precise sense, 
the eigenstates of the initial system will evolve into the 
eigenstates of the final system if this interpolation takes 
place slowly enough. Slow enough in simple quantum 
systems usually means that the variation happens over 
much longer time scales than the characteristic time of 
the system (by which we mean the time associated with 
the typical energy of the system:r ~ i). Having said 
this, were we to study QFT in places where the met- 
ric varies a lot more rapidly (at the edges of black holes 
or in the very early Universe) we invariably have to ac- 
count for the curvature of space. Thus, we can hope that 
the effects of a time dependent background on the closed 
string spectrum only require minor modifications to the 
flat space spectrum, if this time dependence is slow com- 
pared to the characteristic time of the string dynamics. 
We show in the Appendix that this is indeed the case, 
and in what follows we will stay within this regime. 

The Outline of this paper is as follows: we first derive 
the energy-momentum tensor of a string gas (the deriva- 
tion here is more general than the one given in |15(). We 
then insert this tensor into the Einstein equations and 
study the dynamics of the radius of the compact dimen- 
sion, assuming that the three large spatial dimensions are 
in the expanding phase. First, we consider a pure state 
of matter. Next, we extend the discussion to a thermal 
state. In Section IV, we discuss the late time dynam- 
ics and show that the stabilization of the radion is not 
compatible with inflation in the three large spatial di- 
mensions, assuming the simplified description of matter 
which we are using. 

A few words on our notation: Greek indices typically 
stand for 5-dimensional space-time indices, Roman in- 
dices ... are associated with the non-compact spatial 
dimensions, and Roman indices a, 6, ... are string world- 
sheet coordinates. The 5-dimensional Planck mass is de- 
noted by M p i 5 (or M5 in abbreviated form). We also 
work in natural units (c = H = fcs = 1) where we pick 
energy to be measured in electron volts. 

II. THE ENERGY-MOMENTUM TENSOR 

To study how a gas of strings affects space-time dynam- 
ics, we need to derive the energy-momentum tensor of 
such a gas. We begin by studying the energy-momentum 
tensor of a single closed string. Starting with the Nambu- 
Goto action 

S„a = ^f**V=K, (1) 

where h a b denotes the world sheet metric 

h ab = d a X»d b X v g lxu {X) (2) 

(and h is its determinant), we see that any variation in 
the space-time metric induces a variation in the in- 
duced world-sheet metric (where the unmatched indices 
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indicate that we perturb only the A/3 component of the 
metric) : 
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h ab + d a X x d b Xl } 5 D (X T -y T ) 
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Now, varying the Nambu-Goto action with respect to 
the space-time metric (performing a perturbation 5g a p 
which acts on the metric as given above) will give us the 
space-time energy-momentum tensor of a single string: 
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We must first discuss the meaning of the expression 

d 2 a S D (X T -y T ) 

= J d 2 aS(X° - y^X 1 y 1 )..^ - y D ) . 

In order to change the variable of integration, we need 
to apply da a = q X xX and sum over all the zeroes of 

X A [cr] — y x when performing the integration. However, 
since we considering modes winding one particular spa- 
tial direction, there are precisely two coordinates that 
are monotonic functions of the world-sheet parameters: 
X° being a monotonic function of er° and X D being a 
monotonic function of a 1 (the D th direction is taken to 
be compact). Thus, 

d 2 o-6 D (XT[o-]- y T) 

= do-°do-H{X {o-} - y°)5(X l [a) - y 1 )..^ ^} - y D ) 
dX° dX D 
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|X°| \*' . 

xS(X° - y°)S(X D - y D )5 D - 2 {X* - , 

where we include the metric factors in the last line so 
that we can take the delta functions in X° and X D to 
be properly normalized. With this result, we get: 



where we use the inverse metric to write the metric con- 
tributions in the denominator. Thus, the single string 
space-time energy-momentum tensor becomes 
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Inserting the explicit form of the inverse world-sheet met- 
ric 



h ab = ~ 
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and using the constraints on the world-sheet fields 
P^X'» = 

+ X'^X'^ = , 

we can write © as 

-1 S D - 2 {X l -y l ) 
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Next, we solve for X° using the constraint © which 
becomes 
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where we have explicitly used the background metric 

g^ v = diag(-l 7 a 2 {t),a 2 {t),a 2 (t),b 2 (t)). (8) 

It is consistent with the equations of motion in a (slow 
enough) time varying background to set X'° = |31|. so 
that 



X Q = P*P: + X H X' 
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(9) 
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| X° | \*' D \ 

1 
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In addition (in a slowly time dependent background) the 
right hand side can be expressed in terms of the familiar 
oscillator expansion. Accounting for the zero mode op- 
erators explicitly, we get the center of mass momentum 
from the spatial zero modes and the winding energy from 
the zero mode terms in the compactified direction. The 
_ other modesgive us the left and right moving oscillator 

hh ab d a X x d b Xt1rxns (see [13 for details): 
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S u ~ 2 (X l - y l ) 
I X° X' D \ 
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x v = Sjgvpipj + ^(N+N -2) + (~)2 + C^y , (10) 
V a! b a 

where n and w are the quantum numbers for momentum 
and winding in the compact direction, respectively, and 
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N and TV are the levels of the left- and right-moving os- 
cillator modes of the string, respectively. The expression 
above is none other than the energy of the string. Using 
the level matching constraints 



N + nw- N= . 
we finally end up with [8^ | 
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Now, we are ready to evaluate J7J for a single wound 

string. We have an explicit expression for X° and we 
know that \X' 5 \ — \w\b in units of a' for a wound string, 
the factor of \w\ being canceled by the summation over all 
(w in total) zeroes of the argument of the delta function. 
Thus, component by component, we get: 
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(note that we label the extra spatial coordinate by "5" ) 
where we ignore off-diagonal components since we are 
about to apply these expressions to an isotropic gas of 
strings. 

However, we wish to present at this point another 
derivation of this result which is rather more direct. Con- 
sider the energy of a single wound string: 
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A spatially uniform gas of such strings with the same 
quantum numbers would have a 5-dimensional energy 
density 
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where fj,(t) is the number density of strings. We divide by 
2irb since this energy will be uniformly distributed over 
the length of the string. The momentum that appears in 
this expression is now the momentum squared of a gas of 
strings whose momenta have identical magnitudes, but 
whose directions are distributed isotropically. To fully 



account for the metric factors in this expression, we write 
n(t) as fio(t) / a 3 (t) since this is how a number density 
explicitly depends on the metric. Now, realizing that this 
is an energy density, we can introduce this gas of strings 
as matter interacting with the gravitational field by just 
adding the following term to the gravitational part of the 
action: 
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(18) 



(see e.g. Section 10.2 of [20| ) . 

Realizing now that the metric factors in the denom- 
inator of the expression for the energy density can be 
written as a 3 — y/ det(g i j) and b — ^/g^, we can write 
the above equation as: 
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By our metric ansatz and the isotropy of the distribution 

n 2 2 2 

of the momenta, we have that p l pi = a (^- + ^- + ^-). 
Using this fact, it is straightforward to show that the 
energy-momentum tensor derived from this interaction 
term is: 
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(19) 
(20) 
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which is exactly what we would get from (|13fl . (|14|) and 
(I15JI were we to construct a hydrodynamical average with 
an isotropic momentum distribution. 

We now investigate some simple aspects of our result. 
The first thing to note is that T| , which is the pressure 
along the compact direction, gets a negative contribution 
from the winding of our strings and a positive contribu- 
tion from the momentum along this direction. The spa- 
tial pressure is always positive, and for the simple case 
n = W = 0, N = 1, which describes a gas of gravi- 
tons moving in the non compact directions, we obtain 
r = 0,p = p/3. 
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Since we are about to study the effects of this energy- 
momentum tensor on space-time, we should make sure 
that the energy-momentum tensor is covariantly con- 
served, or else it will not be consistent to equate it to 
the covariantly conserved Einstein tensor. The covariant 
conservation of 

= V M T^ , 

where is the covariant derivative operator, implies 



P+3-(p + p) + j{ P + r) 
a b 



dip 
d 5 r. 



It is straightforward to check that our energy-momentum 
tensor satisfies this as an identity. In the continuity equa- 
tion, this is due to the metric factors contained in the 
energy density, which upon differentiation produce terms 
that exactly cancel the terms proportional to the Hub- 
ble factors. The remaining equations are trivially sat- 
isfied by our setup, which assumed an axis of symme- 
try along the compactified dimension (the Kaluza-Klein 
setup) with homogeneous and isotropic spatial sections. 

One final point to note is that we have derived an 
energy-momentum tensor that exhibits positive pressures 
along the non-compact directions and positive or nega- 
tive pressures along the compactified direction. We need 
to ensure that this negative pressure has a bounded equa- 
tion of state as otherwise our theory would be unstable. 
The Dominant Energy Condition (DEC) of General Rel- 
ativity [2l| ensures the stability of the vacuum, and re- 
quires that the equation of state parameter lo = p/ p be 
greater than or equal to -I (see e.g. (2^ for a recent 
discussion). Since the spatial pressures are always posi- 
tive, we only need to check our equation of state for the 
pressure along the compact direction: 
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where the co-efficient of p in the above is our equation 
of state parameter. Were we to consider states described 
by n = ±1, w = —ft, N = 1 (which as we will see further 
on, turn out to be the relevant states that give us sta- 
bilization), this parameter remains bounded as b varies 



-1< w < 1 



(22) 



Thus, we have verified that the spectrum of string 
states satisfies the DEC, and in doing so ensured our- 
selves of sensible space-time dynamics arising from the 
string gas, the topic we will turn our attention to next. 



III. SPACE-TIME DYNAMICS 



We start with the Einstein tensor derived from the 



metric ©: 
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h 
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Equating this to , A T{f will give us the Einstein equa- 



tions. However, let us focus on the equation that governs 
the evolution of the scale factor b. Starting with G* and 

- 2 

eliminating a and a by adding the appropriate combi- 
nations of Gg and G|, we get: 



b+3Hb 



p - 



2r p 

y ~ 3 



0. 



(23) 



where H is the 3-dimensional Hubble factor. This is a 
second order, nonlinear (because of the b dependence in 
the matter terms) differential equation with a damping 
term and a driving term. We will demonstrate further on 
that the Einstein equations admit expanding solutions 
for the non-compact dimensions {H > 0), and take it 
as a given for what follows. Thus, in spite of its non- 
linearity, we easily see that (|23|l describes an expanding 
or a contracting scale factor depending on the sign of the 
driving term. 

The first thing to notice from this equation is that 
matter for which the quantity p — ^ — | vanishes will 
not contribute to the dynamics of the compact dimension. 
Thus, recalling that a gas of gravitons (n = w — 0, N — 
1) has an equation of state p = ^, r = 0, as does a 
gas of ordinary 4-dimensional photons, we see that such 
matter will not affect the dynamics of the scale factor 
b. In fact, such a background 4-dimensional gas provides 
an excellent candidate for a thermal bath which we will 
eventually couple our gas of winding modes to. 

First, however, we will study this driving term as it 
is, for a gas consisting of strings with identical quantum 
numbers. Upon evaluating the driving term, we find that: 



M 3 , 

ph 



P 



2r 



Mo 



(24) 
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4(iV-l) 
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from which we infer that momentum modes and oscillator 
modes lead to expansion of the scale factor, whereas the 
winding modes produce contraction. Exactly what hap- 
pens, of course, depends on the values of the quantum 
numbers. It should be recalled that the quantum num- 
bers are subject to the constraint nw + N > coming 
from the level matching conditions (see Eq. Illfl . 



6 




FIG. 1: Potential term for n = -w = ±1,N = 1. The 
horizontal axis is 6 (in string units), the vertical axis gives 
the potential in units of U, where U is the prefactor on the 
right hand side of H25H . 



Let us pick a particular set of quantum numbers. As 
we shall see later, the most interesting case is when 
n = —uu = ±1, N = 1, in which case the driving term 
becomes: 



At this point, we wish to mention that the "Quantum 
Gravity Effects" required to stabilize the extra dimen- 
sions in earlier attempts [23l l24j at Kaluza-Klein cos- 
mology find a stringy realization here, in that all that 
was required for radius stabilization was matter that de- 
pended on the size of the extra spatial dimensions in a 
non-trivial way. 

To round off the discussion, we wish to demonstrate 
that our assumption of an expanding scale factor a(t) is 
consistent with an oscillating scale factor b{t). Consider 
the two Einstein equations that do not contain second 
derivatives of b: 
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These equations imply that 



II =-^(2p + r). 



(27) 
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'a' + ii~b) 2 
(25) 

where b is the scale factor in units of \fa! . Quite gener- 
ically, we can explore features of the "potential" energy 
that will yield this driving term. We see that, since the 
denominator is strictly positive, and the driving term 
changes sign at b= 1, this value will be a minimum of 
the potential energy, and hence a point of equilibrium. 
Numerical integration of the dri ving term yields the po- 
tential energy curve of Figure 1 j34| where the potential 
is plotted in units of 2 ^°_ j= as a function of b~. 
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Because of the Hubble damping term in the equation of 
motion for b (which is obtained by dividing l|23[) through 
by a factor of \/c7), the scale factor will perform damped 
oscillations about the minimum of the potential to which 
it will evolve with rapidity depending on the value of the 
"spring constant" multiplying the driving term: 



2^0 



M^ h a 3 2Tra' 



(26) 



Thus, we have established that a gas of string modes 
with non zero winding and momentum numbers in the 
compact direction will provide a dynamical stabilization 
mechanism for the radion, provided that the three non- 
compact dimensions are expanding (such a behavior was 
already found in an early study [lOj of the dynamics of 
string winding modes - we thank Scott Watson for draw- 
ing our attention to this paper). We will address the 
phenomenology of this stabilization mechanism further 
on, simply stating for now that we can obtain a robust 
stabilization mechanism which is consistent with obser- 
vational bounds. 



The resulting equation for H has the integrating factor 
1/6 2 , and hence the solution: 



H{t) 



H 



3M| 



(28) 



Now, from the discussion surrounding Ij22(l . we see that 
p > r > —p. Thus, the contribution to the integral in 
the above is strictly positive, and the second term on the 
right hand side of 1|28[) can at most take on the value: 



Ml 



dt 



, Pit') 



(29) 



Thus, we see that if we pick the initial conditions for 
H appropriately, the scale factor a can be taken to be 
expanding (H > 0) regardless of the detailed motion of 
b. In fact, if we assume that Hq starts out positive, then 
H(t) will remain so if 



Hn > 



1 



3A/f 



dt > 2 P + r 



(30) 



where the eventual stabilization of b and the 1/Vol de- 
pendence of p and r will bound the integral, which im- 
plicitly depends on H itself. This implicit dependence 
works in our favor in that the larger we make H a , the 
smaller the integral becomes and so we can imagine pick- 
ing an initial Hq such that a persistent expansion of the 
non-compact dimensions results. Note, however that if 
in the spirit of brane gas cosmology, we assume that all 
spatial dimensions are starting out with the same size 
and instantaneously static, then it may not be possible 
to evolve to a situation in which three large spatial di- 
mensions are expanding. This is, in fact, the result that 
emerges from the work of |13| , at least in a certain region 
of phase space. 
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IV. THERMAL STRING GASES 

In what we have done so far, we have just considered 
the behavior of the size of the extra dimension in a rather 
artificial setting, namely imposing a gas of strings with 
a fixed set of quantum numbers. One expects the early 
Universe to be in a state of thermal agitation, and it is 
inevitable that transitions between different energy levels 
will be induced in the string gas. Thus, to have any hope 
of realistically applying our setup to cosmology, we need 
to study the effects of placing the string gas in a thermal 
bath. Referring to our expression for the energy density 
of a string l|19[) , we see that a gas of strings with different 
quantum numbers will have the energy density: 



p= E 

n,w, N,p 2 



f-t>n,w,N,p 2 

a 3 b2ir 



p* pi+ (JV -1 + - + — 2 

r\:' h rx' 



(31) 

with densities fJ- n ,w N p 2 f° r each given set of quantum 
numbers. The expressions for the pressure terms p and r 
are similarly modified. If we are in thermal equilibrium, 
the densities are given by the Boltzmann weight 



u v 2 = e 0E ^( 



fi re f 



(32) 



where the subscript "ref ' refers to some arbitrary refer- 
ence energy level. 

What constitutes the thermal bath to which the string 
gas is coupled to? We know from the discussion at the 
end of Section II that gravitons described by unwound 
strings propagate in the non-compact directions with an 
equation of state p = p/3. Introducing a gas of ordi- 
nary photons will also add a 4-dimensional component 
to the energy-momentum tensor with the same equation 
of state. Such particles offer us an ideal candidate for 
a thermal bath, for two reasons. Firstly, thermal equi- 
librium demands a coupling of some kind between the 
gas of winding modes and the gas of gravitons and pho- 
tons. Such a coupling is readily provided by the tree- 
level reaction w+ w— > via which winding modes of 
equal and opposite winding scatter to produce 4-d gravi- 
tons. This thermalization mechanism will, at non-zero 
temperatures, create an equilibrium where there will be 
an ever-present non-zero winding mode density due to 
gravitons scattering into winding modes (and vice- versa). 
This thermal bath has the further property that it does 
not affect the dynamics of the extra dimension other than 
through the Hubble factor (which it influences), since 
the driving term is only sensitive to the combination 
p — p/3 — 2r/3 which vanishes for the graviton and photon 
components of the energy-momentum tensor. 

With the above in mind, the driving term for the scale 
factor b becomes: 



/ 2r p 

l p -y-3 



e /SE rcf 



(33) 



E 

n,w,N,p 2 



-/3E 



b 2 



2nw 
3a' 



3a' 2 



4(iV-l) 
3a' 



VVi 



-!) + ($ + a*) 



b\2 



where the Boltzmann weight in the summation depends 
on the values of the quantum numbers. We remind the 
reader that the sum is restricted by the level matching 
condition N + nw > 0. For completeness we also remind 
the reader of the resulting equation of motion for b: 



= b+3Hb 



PE ref 



M pl 5 a32lT 



(34) 



E 



-PE 



,N,p 2 



f n 2 



2nw w 2 b 2 4(iV - 1) 



3a' 



3a' ; 



3a' 



The summation which has to be performed in order 
to obtain the driving term is quite formidable, were it 
not for a rather special feature of string thermodynam- 
ics. Consider the argument of the exponential in the 
Boltzmann factor: 



4 , K ,n wb,„ 
0E n ,w,N,p 2 = (3\/p i Pl + — (N- !) + (- + — ) 2 



(3 



a'p l p t + A{N - 1) + (-=■+ b w) 2 



We see that when the energy is expressed in terms of di- 
mensionless variables, we pull out a factor of \fot ' . Thus, 
the argument of the exponential in the Boltzmann weight 
is P/y/a' times a term of order unity. To be able to ne- 
glect all but the first few terms in the summation, we 
need the Boltzmann factor to be considerably less than 
unity, i.e. that 



< 1 



Thus, if this condition is satisfied, then the terms which 
dominate the sum will be those whose quantum num- 
bers render them nearly massless, since any state with 
even one of its quantum numbers being different from 
the nearly massless combination will produce a term of 
order unity times (3/\fa'. 

Let us take a closer look at the above condition. We 
know from string thermodynamics that there exists a lim- 
iting temperature - the Hagedorn temperature |25| . 
Thus, for us to even be able to apply thermodynamics, 
we need to be well below this temperature, which for all 
the string theories is of the order of Th ~ l/\fa' \ Thus, 
~ vc?, and so if we are at a temperature much lower 
than the Hagedorn temperature, i.e. T <C Tjj or equiva- 
lently (3h <§C (3, then 



r a l < f3. 



(35) 



which is exactly what we need for the Boltzmann weights 

of higher mass states to be negligible. So, even if the 

thermal bath has a temperature of only one order of mag- 

§_ 

nitude below the Hagedorn temperature, then e vST 1 ~ 
10~ 5 which clearly lets us ignore any term whose energy 

in dimensionless units ^Ja'p l Pi + 4(iV — 1) + (j+ b w) 2 

is anything other than zero. This translates into us being 
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able to neglect all states other than those that are mass- 
less. Thus, the summation now becomes very tractable, 
and we can also have faith in our truncation of the string 
spectrum to the lightest states all the way up to very 
high temperatures (T ~ 7^/10). Before we carry on we 
should remark that exactly massless states have a non- 
zero momentum given by the thermal expectation value 
of E = \p\ = 3/f3. 

Let us then proceed to evaluate l|3*3*|> . so that we can 
evolve b in time using (|34[) . recalling that now we only 
sum over the massless and near massless states subject 
to the level matching constraint. Let us begin near b= 1, 
i.e. b= 1 + r. Then for the case that L ^ 1, we only have 
one truly massless state: n — w — 0, N = 1. This term 
will not contribute to the driving force for b since 



2nw 
3a' 



3a' : 



4(Af-l) 
3a' 



\P\ 



= 0. 



(36) 



Thus, the next lightest state which has quantum numbers 
N = l,n = — w = ±1 will dominate the evolution of b. 
The level matching constraints N + nw > ensure that 
there are no more nearly massless states (Note we only 
consider states with positive mass squared - any tachy- 
onic states are posited to be absent from our spectrum). 
Such states will contribute: 



N^ 2 



-GL p 



1 

b 



a- b) 2 + oi'v 2 



(37) 



Expanding b as 1 + L and ignoring terms higher than 
quadratic in T results in a contribution to Ij37(l of: 

Since there are two such terms which yield identical 
contributions, the sum total of the contributions from 
the near massless states yields the equation of motion 



L +3H r 



M 5 3 27ra 3 |p|a' 3 / 2 



(38) 



where the exponential factor gets cancelled if we use this 
massless state as our reference state, as in (|32() . The form 
of this equation clearly shows that T will tend to zero if 
it starts out on either side of this value. 

However, to confirm that L = is a genuine point of 
equilibrium, we need to confirm that the extra massless 
states that appear at this radius (8 in all) contribute 
in such a way so that their sum vanishes. This can be 
verified by a straightforward calculation [35| . 

However, we wish to point out that as long as we stay 
in thermal equilibrium with the graviton gas, this equi- 
librium is actually metastable. The reason for this is easy 
to see from the formula for the mass of a winding mode: 



a m 



^+wb) +4(iV-l) 



(39) 



In addition to the massless state given by n = w = 0,N = 
1 (the graviton), and the 8 other massless states that 
appear at the self-dual radius (which are given by N — 
l,ra= -w = ±1; N = 0,w = n = ±1; N = 0,w = 0,n = 
±2 and N — 0,w — ±2,n = 0), there are additional 
massless states at further special radii 



b= -. — r ; w = ±m, n = 0, N = 
m 



b= ~y ; n = ±m i w = 0, N = 



m e Z . 



Thus, at half-integer multiples and and half integer frac- 
tions of the self-dual radius, two massless modes appear 
and will thus yield the dominant contributions to the 
driving term. These contributions again exactly cancel 
at twice the self dual radius, and at half the self dual ra- 
dius. However in general, we will get a driving term that 
yields expansion at half integer points above twice the 
self dual radius and similarly, contraction at half integer 
fractions below half the self-dual radius. However since 
we posit that we begin at or near the self-dual radius, 
we are guaranteed to stay locked near it if our initial 
conditions satisfy 




M 3 a 3 27ra' 



where the last condition constrains the initial "veloc- 
ity" of the scale factor to be such that it cannot roll 
over the "hump" in the potential energy surrounding the 
metastable equilibrium at b ~ ycV . 

Thus, we have demonstrated in the context of GR 
how a string gas in thermal equilibrium with a bath 
of gravitons and photons will dynamically stabilize the 
scale factor of the compact direction if we begin close 
to that radius. Thermal equilibrium with the graviton 
bath ensures a persistent non-zero density of such wind- 
ing modes. One can now imagine that, at some point, the 
winding mode gas becomes decoupled from the graviton 
gas, i.e. falls out of thermal equilibrium. In this situa- 
tion, we are left with an unchanging driving term of the 
form H24J) . which yields the potential in Fig. 1, which 
will guarantee radion stabilization at the self dual radius 
for the remainder of the Universe's dynamics. We now 
turn our attention to the possible connection between 
this mechanism and inflationary and standard Big Bang 
cosmology. 



LATE TIME EVOLUTION B. Matter Dominated Evolution 



Recall the Einstein tensor for our metric setup: 

Go = + 

a la bi 

J J I a b \a/ o a\ 

r rd /a\ 2 l 

Gt = -3 - + (- • 
La Va/ J 

We know that the dynamics of the scale factor b in the 
situations we studied above cause it to undergo damped 
oscillations around the self dual radius. We demonstrated 
in a previous section how the "spring constant" of this 
evolution will lock in to this equilibrium fairly rapidly. 

We can then study the evolution of the Universe after 
radius stabilization, which implies that b=b= and p — 
2r/3 — p/3 = 0. The resulting Einstein equations are: 



^10 1 rriO 


+ H 2 


P 
3 




( '' uT* 


-> H= 


■> 


hp) 




-> P~ 


y ~ 2 


= 


= - 


■+ P4- 


3ff(p + 


p) = o, 



where the 55 equation is precisely the equilibrium condi- 
tion. Thus, we recognize in the above the basic equations 
of FRW cosmology. We now consider how one achieves 
the two important epochs of late time FRW cosmology, 
namely the radiation dominated era and the matter dom- 
inated era. 



A. Radiation Dominated Evolution 

If we assume that the density of 4-d matter gas is far 
greater than the density of the winding mode gas, then 
this will be the dominant component that drives the evo- 
lution of the macroscopic dimensions. If the 4-d matter 
has an equation of state parameter w, then the solutions 
to the Einstein equations become: 

p{t) oc a-^ 1+w) (40) 

or 

a(t) oc ^ . (41) 

Thus, for a 4-d graviton and photon background, wc 
get that a(t) oc t 1 ' 2 , and so we reproduce a late time 
FRW evolution that is consistent with standard Big-Bang 
cosmology immediately after the end of inflation, whilst 
maintaining radius stabilization. 



Reconsidering J33J): 

b ( 2r p\ 

We see that any matter with the equation of state of 
non-relativistic dust (p = 0), can only drive the expan- 
sion of the radion if it is of a 3-dimensional nature (i.e. 
r = 0). This is surely to be a cause for concern when 
considering that at late times, one (naively) expects 3- 
dimensional non-relativistic dust to be the major driving 
component of the expansion of the universe, which would 
normally invalidate our stabilization mechanism in the 
present epoch. 

However we wish to remind the reader that present 
day observations demand that a significant fraction of 
the energy density of the universe, which also drives the 
present day matter dominated expansion, be in the form 
of cold dark matter- whose nature is as of yet completely 
unknown. There is a significant amount of interest is the 
prospect that extra-dimensional matter or extra dimen- 
sional effects might account for the 'missing' matter in 
the universe. In what follows, we find that the only way 
to make our stabilization mechanism consistent with a 
'matter dominated' epoch is to introduce extra dimen- 
sional cold dark matter. We propose a candidate for 
this cold dark matter which is naturally contained in our 
framework, and discuss other possibilities which might 
have a natural realisation within the general brane gas 
framework (note that a similar proposal was made in 
0). 

We see that in order for the driving term in l|23|) to 
correspond to a stable minimum at the self dual radius 
for matter which obeys the equation of state for non- 
relativistic dust (p = 0), we need to consider matter for 
which: 

r = -p/A = ^- (42) 

That is, we require the dominant component of the en- 
ergy density which is driving the expansion of the uni- 
verse be such that it preserves the stability of the ra- 
dion at the self dual radius. Matter which exhibits such 
an equation of state will surely have to be massive (else 
there will be a non-zero pressure along the non-compact 
directions for any non-zero energy density). In addition, 
such matter will have to be something beyond presently 
supposed dark matter candidates (WIMPS, supersym- 
metric relics etc.) as it will neccesarily have to be 'extra- 
dimensional' in nature. We now show that our model 
naturally contains such a candidate. Recalling the dis- 
cussion surrounding l|22|l . we see that the equation of 
state parameter for a gas of strings with a particular set 
of quantum numbers is obtained from the following equa- 
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tion: 



P x 



&2 



pip i + A T (N-l) + {^ + ^) 



b\2 



(43) 



where the momentum will be set to zero (or is vanishingly 
small) in order to satisfy p = (c.f. 120|l ). In particu- 
lar, since we are looking for states which can satisfy l)42|) 
at the self dual radius, we need to find the appropriate 
quantum numbers which have an equation of state pa- 
rameter w = — 1/2 at 6 = VcV, which reduces to the 
following condition: 



3n 2 - w 2 + i(N - 1) + 2nw = , 



(44) 



and we have to be mindful of the level matching con- 
straint: N + nw > 0. As expected, it turns out that the 
massless states that satisfy these conditions have an en- 
ergy density proportional to \p\, whereas the pressure is 
proportional to |p| 2 /3, and hence one cannot have a non- 
zero pressure without having a vanishing energy density. 

The first massive states which satisfy l|44[) are repre- 
sented by the quantum numbers: 



\p\ = , N = 2 ; n = , w = ±2 . 



(45) 



These states contribute to the energy momentum tensor 
as follows (c.f. (O - ©): 



P = , p = 



2fiy/2 
2na 3 a' 



-p/2, 



(46) 



where the factor of a' in the denominator comes from 
two factors of \fcH- one from the metric factor b which 
has stabilized at the self dual radius, and the other as 
the pre-factor of the non zero rest mass of this string 
state. As we will see in the next section, were we to look 
at fluctuations of the radion around the self dual radius, 
these states also provide a stable equilibrium at b = \foJ 
in a phenomenologically acceptable manner. Thus, tak- 
ing questions of consistency with observation for granted 
for the moment, we see if these particular string states 
are taken to dominate the present energy density of the 
universe, then by the Einstein equations derived at the 
start of this section, we can admit an epoch of dust driven 
FRW expansion (p — 0, a cx i 2//3 ) at late times, consistent 
with radius stabilization. 

However, there are many further issues that will have 
to be resolved if we are to take this idea of stringy dark 
matter seriously, which we postpone to a future report. 
At present, however, we wish to state that there are 
indications that such stringy dark matter might have 
the right clustering properties at the level of first order 
perturbation theory, in that local overdensities of this 
stringy dark matter induces gravitational clustering in 
the remaining 3-dimensions. 

We also wish to point out that in certain situations, 
non-relativistic p-branes might also be able to provide 
us with a matter content with satisfies (|42f) jl5j . In the 



context of brane gas cosmology, this is an appealing idea 
as one might need higher dimensional branes to stabi- 
lize compact sub-manifolds that do not admit topological 
one-cycles (and hence do not admit winding modes). We 
will investigate this possibility further in a future report. 
Finally, we wish to address the effects of an intermediary 
epoch of scalar field driven inflation on our stabilization 
mechanism. 



C. Intermediate (Non-Stringy) Inflation 

We find that our mechanism for radius stabilization 
might not be compatible with an intermediate epoch of 
bulk scalar field driven inflation. To investigate this, we 
first consider the energy-momentum tensor of an almost 
homogeneous inflaton field: 

.2 

TS = Mag (- [t. + vfr)] , [t. _ V{(j)) ] , - V (</>)] ) 

(47) 

Adding this energy-momentum tensor to our string gas 
yields a non-trivial contribution to the driving term in 
the equation of motion for b. This driving term then 
takes the form 



e PE ret 



3M, 



n,w,N,p 2 



Ml a 3 2n 



-HE , „2 2nw w 2 b 2 



n 



3a' 3a' 2 



(48) 

4(N-iy 

3a^ 



from which it is easy to see that the inflaton contribution 
will drive expansion in the extra dimension. In general, 
this term will compete with the string gas contribution 
which, as we have seen, drives contraction, if we are above 
the self-dual radius. However, this competition is short 
lived, as the factor of a 3 in the denominator of the string 
gas driving term will quickly render it irrelevant and the 
scale factor will then expand according to 



b+3Hb 



2b 

3M 3 



V{ct>) = 0. 



(49) 



Recalling that during this (slow roll) inflation H, and by 
the time-time Einstein equation, also V{4>) are almost 
constant, we can solve the above equation, with the re- 
sulting two solutions: 



bit) (x exp 
b(t) oc exp 



,l + ,/l + 7 



Substituting in the Einstein equation H 2 = V((f>)/3M§ 
gives us 



b{t) oc exp-V 
b(t) oc exp~ 



oc e 
oce 



-3.56Ht 



I v I ■:-§')* , n.-i(i/// 
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Except for very special initial conditions, the growing 
mode will rapidly come to dominate. Thus, we conclude 
that b expands exponentially (though not as fast as a). 
After inflation has finished (and after reheating to a tem- 
perature smaller than the one required for the pair pro- 
duction of string winding modes), a expands as t 1 / 2 . The 
energy density in the string gas will have been exponen- 
tially suppressed by the inflationary evolution, and thus 
to good approximation the equation of motion for b will 
take the form 



leading to 



which implies 



b +3H b= 



b= Ca- 3 



b(t) =6(1) + 2 6(1) 1 



(50) 



(51) 



where (1 < t). Thus, b asymptotically expands to some 
limiting, and very large (due to the initial conditions that 
result at the end of inflation) value. 

In conclusion, we have seen how our radion stabiliza- 
tion mechanism is consistent with the FRW evolution of 
the non-compact dimensions, but not with an interme- 
diate inflationary period, with inflation driven by a bulk 
scalar matter field. Thus, in order for brane gas cosmol- 
ogy to make contact with the present cosmological ob- 
servations, one either needs a different (maybe stringy) 
realization of inflation (see f° r some ideas) where 
strings are produced in re-heating, or a non-inflationary 
mechanism to solve the cosmological problems of stan- 
dard Big Bang cosmology and to produce a spectrum 
of almost scale-invariant cosmological perturbations. Fi- 
nally, we turn to various phenomenological issued per- 
taining to this model. 



VI. PHENOMENOLOGY 

There are two potential phases of applicability of our 
considerations. The first is to the early phase of string gas 
cosmology before a period of inflation. In this case, there 
are no phenomenological constraints on the model since 
the number density of the particles (from the four dimen- 
sional perspective) which correspond to the string states 
wrapping the extra dimension are diluted exponentially 
during inflation. However, in this case our considerations 
would no longer be relevant for the late-time stabilization 
of the extra dimension. 

The second phase of potential applicability of our con- 
siderations is to the universe after inflation of our three 
spatial dimensions. We then need to assume that winding 
and momentum modes about the extra dimension can be 
regenerated in sufficient number, as discussed in pr| . In 
this case, there are important constraints on our model. 



We must ensure that the particles corresponding to our 
string states do not overdose the universe. In addition, 
there is a radion mass constraint. Since the radion cor- 
responds to a scalar particle from our four-dimensional 
perspective, we must make sure that its mass is consistent 
with the experimental constraints (we thank the Referee 
for stressing this point to us). 

An additional constraint comes from the string the- 
oretical aspect of our model: we must ensure that the 
derivatives of the metric remain several orders of mag- 
nitude smaller than the worldsheet derivatives (see Ap- 
pendix). This is to ensure that we can inherit aspects of 
the string spectrum and constraint algebra that we have 
used so far. From l|38|) . which has the form: 



T + ZHT + kT = 



(52) 



we see that the 'spring constant' which sets the scale for 
the how fast the metric factor b varies, is given by 



k 



8^o 



3Mla' 3 / 2 \p\a 3 



(53) 



where the subscript on /i is to remind the reader that any 
explicit metric dependence has already been factored out 
(see section II). We have taken the stabilization to be 
provided by the massless states discussed earlier with |p| 
denoting the momentum in the non compact directions, 
N = 1, and n = — w = ±1. In order that our metric fac- 
tors evolve much slower than the string scale, we require 
that d^g <C 1/ \j2ixa! '. Since (jSHJ is a second order ODE, 
this implies that k <C l/2wa'. As discussed in the Ap- 
pendix, we choose to be quite conservative and demand 
that k < lCT 6 /27ra'. 

A second constraint comes from requiring that the sta- 
bilization mechanism be effective at all times. This leads 
to a lower bound on k. We take this bound to be given 
by the 'critically damped' value for k: 



m 1 



(54) 



The two above constraints yield the following bounds: 



9H 2 < 



3M 5 3 a' 3 / 2 |p|( 



< 



KT 



27ra' 



(55) 



Since the winding states that stabilize the extra di- 
mensions are massless at the self-dual radius, they will 
behave as hot dark matter - dark matter because they 
only interact gravitationally (through the tree level in- 
teractions w+ w— > h^v) with other fields, hot because 
they are massless and have a radiative equation of state. 
We have to ensure that we do not introduce too many 
of these objects so that we can ensure consistency with 
observational bounds. 

Next, we have to ensure that the massive string states 
that we propose as a candidate for the cold dark matter 
that is presently driving the 'dust dominated' expansion 
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of the universe do not violate any observational bounds 
while preserving our stabilization mechanism. 

Let us begin by considering the massless states which 
are presently stabilizing the radion. From 117f) we see 
that post stabilization, for stringy matter with N = 1, 
n = — w = ±1, we must have: 



Molgl ^ in -4 

P = 5— < 10 Pcrit 



(56) 



in order to ensure consistency with the nucleosynthesis 
bounds. The critical density of the universe is p cr it = 
10 _29 <//cra 3 . We then find that |5^)l becomes 

Ho < KrtCT^GeVVr 1 ■ (57) 

Let us now parametrize the present value of \p\ as 

|p| = l(T 7 eV (58) 

where 7 is some constant determined by the initial 
conditions and the history of the universe, and this 
parametrization being motivated by the fact that \p\ is 
likely to be of the order of a few eV in the present epoch 
if it corresponds to an initial \p\ of the order of the Planck 
energy. Then, the above bound takes the form 

Po < l(T 41 GeV 3 10 7 . (59) 

On the other hand, the first inequality in l|55|) becomes: 



Mo 



> H 2 {t )\ P \ ~ KT 93 l(T 7 GeV 3 



(60) 



which is consistent with (|59|l . In the above, we have made 
use of 



M 5 " 3 



8ttG 5 



8ttGa 



(61) 



Using H61fl . it can easily be checked that the upper bound 
on po which follows from the second inequality in l|55(l is 
much weaker than the bound H59|) . We conclude that one 
can easily arrange the number of string modes such that 
stabilization of the extra dimension is ensured and at the 
same time the massless modes do not conflict with the 
nucleosynthesis bound. 

Let us next turn to the radion mass constraint: Since 
the radion appears in four dimensions as a scalar field, 
its mass must be larger than 



m crit = 10~ 12 GeV 



(62) 



in order to avoid fifth force type constraints. Since the 
square of the radion mass is given by fc, this constraint 
becomes 



Po > Mloi 



3 n '3/2 



\p\i 



10 _ 35 - 7GeF 3 



(63) 



which is consistent with the upper bound (|59|l on pq if 
7 > 3. Such a value of 7 is not at all unreasonable and 
could easily arise from an additional suppression of the 
momentum during a short period of inflation. 



It turns out to be crucial that we use massless modes 
to stabilize the extra dimensions, as more massive string 
states would bring down the upper bound, and as we are 
about to see, do not provide as effective a spring constant 
and hence bring up the lower bound, to the net effect that 
it is phenomenologically inconsistent to have these as the 
only strings that are stabilizing the radion. We arrive 
at this observation by considering the second aspect of 
our phenomenology- namely, that we would like the cold 
dark matter content of our present universe to consist of 
massive string modes (which satisfy (|42f) . which as we 
have seen is require in order to maintain stabilization at 
the self dual radius). 

Considering the contribution to the energy density by a 
string gas with quantum numbers \p\ = 0, N = 2, n = 0, 
w = ±2 (|46f) . and equating this to the critical density of 
the universe, we see that: 



2y/2p 



dm 



2 7 ra 3 a' 1 /2 



Pcrit 



(64) 



Where now the subscript on p serves to indicate that this 
is our dark matter candidate. This requires 



Pdr, 



10- b7 GcV 3 



(65) 



Upon perturbing around a stabilized radius, we find that 
these string modes contribute to the stability of the ra- 
dion with the spring constant: 



^dm 



ft dm V8 

27ra'a 3 3 



(66) 



Demanding that this value of k is consistent with the 
radion mass bound yield a lower bound on pdm 



Pdr, 



,2 

l cnt ± 



10" 6 GeV 



(67) 



which is clearly inconsistent with l|64|) for values of 7 
which are extremely large. 

Thus we see that if we introduce the correct amount 
of our dark matter candidate, it contributes too feebly to 
the dynamics of the radion (even though it does provide 
its own contribution to stabilization) . However, this is of 
no concern to us, as we have already shown that the mass- 
less string states provide a robust stabilization mecha- 
nism that is consistent with observational bounds. Thus, 
if our string gas has a massive component that serves as 
cold dark matter, and a massless component that sta- 
bilizes the radion (and behaves like hot dark matter) in 
the right proportions, which as we have shown is quite 
easy to achieve, we can be assured of the phenomenolog- 
ical consistency of our stabilization mechanism with late 
time FRW cosmology. 



VII. CONCLUSIONS 



The analysis in this paper is motivated by brane gas 
cosmology [HQ- As a simplified problem, we have stud- 
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ied the effects of a gas of strings with non- vanishing mo- 
mentum and winding modes about a single compact ex- 
tra dimension taken to be a circle (the three large di- 
mensions are taken to be spatially flat and isotropic) on 
the evolution of the radius of that dimension, assuming 
that the background space-time satisfies the equations of 
motion of General Relativity. We discovered that such a 
string gas leads to a dynamical stabilization mechanism 
for the radius of this dimension, the radion. Assuming 
initial conditions in which the three large dimensions are 
expanding, we found that the radion performs damped 
oscillations about the self-dual radius. 

In a first step, we studied the effects of a gas of non- 
interacting strings, each string having the identical mo- 
mentum and winding numbers. Key to the stabilization 
mechanism is the fact that winding modes and momen- 
tum modes contribute with opposite sign to the pressure 
of the string gas in the direction of the compact dimen- 
sion, and that the winding modes generate a potential 
for the radion which favors contraction, whereas the mo- 
mentum modes generate a potential favoring expansion. 
We then showed that the stabilization mechanism also 
holds for a gas of strings in thermal equilibrium. 

We also showed that, after radion stabilization, the 
scale factor for the three large spatial dimensions obeys 
the usual FRW equations of standard Big Bang cosmol- 
ogy. Thus, our scenario leads in a natural way to a 
late time FRW Universe. However, we have also shown 
that the radion stabilization mechanism is not compat- 
ible with a period of scalar field driven bulk inflation. 
Thus, in order for brane gas cosmology to make success- 
ful contact with present cosmological observations, one 
cither needs to find a stringy mechanism for driving infla- 
tion where strings are produced in the re-heating epoch, 
or else one must provide an alternative to inflationary 
cosmology, both for solving the mysteries of standard Big 
Bang cosmology, and for explaining the origin of the ob- 
served large-scale fluctuations. 

Note that we start with the assumption that three spa- 
tial dimensions are already much larger than the other 
ones (one other dimension in our case). Whether or not 
the dynamics of strings in the initial stages will indeed 
lead to this situation may depend on the corner of M- 
theory one is working in, i.e. on the specific form of the 
background equations of motion and initial conditions 
(see 0, 0, H, 0, El for different angles on this issue). 
However, it should not be hard to generalize our analysis 
to a situation with more compact dimensions or different 
topologies and geometries of the extra dimensions, which 
will be the focus of our future work. 

In future work we also plan to study the annihilation 
rate of the string modes which are central to this work, 
namely modes which have both winding and momentum 
in the compact direction. Since these modes interact 
only gravitationally just like gravitons, they will be out 
of thermal equilibrium at late times and hence will not 
decay. 
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APPENDIX A: THE STRING SPECTRUM IN A 
TIME DEPENDENT BACKGROUND 

Let us then begin with the Polyakov action for a closed 
string: 

S= 4^f ^V^^X^X^pO . (Al) 

Varying this action with respect to the world sheet metric 
gives us the equation of motion 

lab = Xh ab = \d a X»d b X»g^ (X) , ( A2) 

where we can exploit the world-sheet diffeomorphism and 
Weyl invariance to make the world-sheet metric flat (con- 
formal gauge): 

la b = l ab = ( "q 1 j) • (A3) 

Varying the action with respect to the world-sheet fields, 
and imposing the gauge choice yields the equations 

d a (d a X x g Xt2 (X)) = ^d a X x d a X v g Xv ^{X) , (A4) 

where the meaning of the derivatives of the metric should 
be clear. This equation translates into 

d a d a X T + T T Xl ,d a X v d a X x = . (A5) 

Now, we consider the case when our metric depends 
only on time, that is g flu (X) = g flu (X°), is diagonal and 
has -1 as its 00 entry (This last point is not essential to 
the argument, it only serves to simplify the equations). 
Unpacking the above equation yields the equations 

-d 2 T X° + d 2 a x° = -^ gxufi (-d T x»d T x x + d a x v d a x x ) , 

(A6) 

and 

-dfrX* + d 2 a X l = -g u g ufi (-d T X l d T X° + d a X l d a X°) . 

(A7) 

On the right hand sides of the equations there is an over- 
all factor containing the time derivative of the metric. 
To estimate the magnitude of either side of the equa- 
tions, we realize that world-sheet time derivatives will be 
of the order of the energy of the string, which is of the 
order of the square root of the string tension: d T ~ — 3= . 
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Similarly, the world-sheet spatial derivative will be of the 
order of the inverse of the string length l s , d a ~ j-, which 
is of the order of the square root of the string tension. 
On the right hand sides, the terms are of the same order 
as the terms on the left, except for the factors contain- 
ing derivatives of the metric. As long as these metric 
derivatives are several orders of magnitude smaller than 
the string tension scale, we can safely ignore them |36| . 
Assuming such a background, l|A6|l and i|A7fl reduce to 
the usual flat space wave equations, and we can proceed 
to expand the solutions in terms of plane waves. We 
then impose the canonical commutation relations on the 
expansion coefficients [37) . 

To complete the analysis, we impose the constraints 
coming from our gauge choice: 

Kb = ^ J) (A8) 

d T X»d T X»g^(X)+d„X>*d„X»g^(X) = 0,(A9) 
d T X»d<,X v gilv {X) = 0(.A10) 

By implementing our plane wave expansion for X^ in 
terms of the creation and annihilation operators, we ob- 
tain the spectrum by requiring that physical states are 
annihilated by the constraint operators. However, when 
we write the constraints in terms of the canonical vari- 
ables, their metric dependence becomes manifest: 

P^X'v- = (All) 



P^P* + X'^X'^ = . (A12) 

We see that (IA11|I is independent of the metric. Since 
this constraint is the origin of the Virasoro algebra, we 
see that it remains valid in suitably time-dependent back- 
grounds. However l)A12jl does depend on the metric since 
P M requires the inverse metric to raise its index and X'^ 



requires the metric to lower its index. Now we choose to 
work in a 5-d background with metric: 

g»u = diag(-l, a 2 {t), a 2 (t), a 2 {t), b 2 {t)) , (A13) 

where the 5'th dimension is taken to be compact with ra- 
dius 2nb With this as our background metric, (|A12|I 
becomes (for a string wound along the 5'th dimension): 

/ - • ;/'•>,/>., + ^{N+N-2) (A14) 
+ <r rrJ'. I Y'\Y' ;; . 

where all we have done is expanded out ilA12(l . and re- 
alized that the terms coming from the non-compact X 11 
and P^ give us the center of mass momenta and the left 
and right oscillator pieces, and the terms coming from 
X 5 have been accounted for explicitly. We know that 
this part of the energy contributes [H| : 

so that combined with the level matching conditions 

nw + N- N= (A16) 

we get 

E = \j 9 ^P l P ] + ^{N-l) + C{ + ^) 2 , (A17) 
V a b a 

where the only remnant of the level matching condition 
is the requirement that 

N + nw > . (A18) 

Thus, we wee that the only effect of working in a slowly 
varying background is to introduce time-dependent met- 
ric factors in the obvious places in (|A17(I which is other- 
wise identical to the result we would obtain in a static 
background. 
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Examples are particle creation, non-uniqueness of the 
vacuum, non-trivial issues concerning existence of asymp- 
totic states. 

Which come from working in conformal gauge h a t = 
diag(—l, 1) 

See the Appendix for all statements made in this sec- 
tion concerning results that are valid in a time dependent 
background 

See the second footnote in the Appendix where we remind 
the reader why the momenta must be contracted with the 
inverse metric. 

Where key to this is the observation that as we approach 
the value b = VcV, these states become massless, and 
acquire a non-zero momentum along the non-compact 
directions, which depends on the ambient temperature. 



If this momentum is large enough, we can be assured that 
I22H is satisfied. 

We model the momentum squared as a smooth function 
of b such that it takes on some non zero value at b= 1 and 
falls off on either side. This is because at b= 1 the state 
described by n = — w = 1, N — 1 becomes massless and 
should have a finite non zero momentum, but as the scale 
factor increases or decreases the state becomes more and 
more massive and hence the momentum becomes negli- 
gible in comparison. The generic feature of a minimum 
at 6= 1 is robust, however, since as we mentioned earlier, 
the force term will always change sign at b— 1. One finds 
that the minimum is always fairly concave independent 
of the nature of the b dependence of the momentum. 
The only non-zero contribution to the driving term comes 
from the states n = 0, w = ±2, N = and n = ±2, w = 
0, N = which make equal and opposite contributions 
and hence cancel. 

For instance, we are tempted to be conservative and 
to ask for them to be roughly 6-8 orders of magnitude 
smaller, in order to be certain of a consistent treatment. 
Recall that d T X^ := P M is the canonical momentum, 
where we must view P M as having its index lowered which 
is to be raised with the inverse metric. Not realizing this 
will produce nonsensical results elsewhere (such as in 
the computation of the energy-momentum tensor of the 
string) in addition to making quantization very awkward. 
This fact is easier to understand if we recall that P M and 
X 1 * are canonically conjugate world sheet fields and not 4- 
vectors. With this in mind, canonical quantization means 
imposing [P m (t, <t),X u (t, <t')] — —iS^8(o — a') which, as 
is, does not involve the metric. Hence we do not have to 
do anything different at this stage of the analysis. 
The results in our paper easily generalize to backgrounds 
of any number of non-compact dimensions so long as pre- 
cisely one dimension is compactified on a circle. For con- 
sistency, we could state that our background is a com- 
pactified 10 dimensional space with six dimensions com- 
pactified on a Calabi-Yau (CY)space and one dimension 
compactified on a circle. Since in the prototypical com- 
pactification scenarios such as the Horava- Witten model 
|27j| the radius of the CY is smaller than the radius of 
the circle, since we will be interested in string winding 
modes but Calabi-Yau spaces do not admit one-cycles, 
we can ignore the presence of the Calabi-Yau space if 
we work in an effective Lagrangian description valid at 
energies smaller than the energy scale of CY compactifi- 
cation. We could always go back and work in a 9 + 1 x S 1 
or 24+ 1 x 5* 1 space-time where we would derive the same 
conclusions as we do here. 



